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^ ABSTRACT 

We examine Gaussian isocurvature hot dark matter (massive neutrino) 
models for large-scale structure in which the initial density perturbations are 
^ produced in the baryons with a power-law spectrum Pb(/c) = Ak" B . We 

^ calculate the linearly-evolved power spectrum and cosmic microwave fluctu- 

ations. We find that models with only isocurvature perturbations are in- 
consistent with observations of damped Lya systems and COBE constraints 
on the power index. However, models which contain a mixture of adiabatic 
and isocurvature perturbations can be made consistent with COBE, galaxy 
surveys and damped Ly« systems. Isocurvature hot dark matter models also 
produce a bias between baryons and neutrinos even in the linear regime. We 
find that this "natural bias" can increase the baryon fraction in small scale 
objects like damped Lya systems, but it has no effect on cluster scales. 
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1 Introduction 



One candidate for the dominant matter in the universe is the massive neu- 
trino, with a mass in the eV range. Since a neutrino with a mass in this range 
is still relativistic when structure formation begins at the epoch of equal mat- 
ter and radiation, massive neutrinos act as hot dark matter (HDM) for the 
purposes of structure formation. Adiabatic (curvature) hot dark matter mod- 
els with a Harrison-Zel'dovich power spectrum suffer from a well-known lack 
of power on small scales, due to the free-streaming of the neutrinos (White, 
Frenk, and Davis 1983; Kaiser 1983; White, Davis, and Frenk 1984; for an 
opposing point of view, see Melott 1985). This problem can be ameliorated 
by beginning with isocurvature fluctuations, either Gaussian fluctuations in 
the baryon component (the subject of this paper), or non-Gaussian isocur- 
vature fluctuations due to some relic "seed" object. In these models, the 
isocurvature fluctuations preserve small scale perturbations until the neu- 
trino streaming length becomes small enough to allow perturbation growth 
on such small scales, significantly enhancing the power on these scales. 

The non-Gaussian HDM seed models have been investigated in great de- 
tail in recent years. Models with HDM and cosmic strings (Bertschinger 
& Watts 1988; Scherrer, Melott, & Bertschinger 1989; Albrecht & Stebbins 
1992) or HDM and generic "seeds" (Villumsen, Scherrer & Bertschinger 1991; 
Gratsias et al. 1993) can produce much more power on small scales than adi- 
abatic HDM and can give much better agreement with the observations. Hot 
dark matter with global texture (Cen et al. 1991) more closely resembles adi- 
abatic HDM, because the texture unwinds on a relatively short time scale 
and does not remain around to seed structure at late times. 

Here we examine Gaussian isocurvature HDM (IHDM) models, in which 
the fluctuations are generated initially in the baryon component. Gaussian 
models have the advantage of simplicity; the initial power spectrum totally 
specifies the model. Gaussian isocurvature perturbations are not as well- 
motivated as adiabatic perturbations, which arise naturally in the context of 
inflation. However, a number of models have been proposed which can give 
rise to isocurvature fluctuations in the baryons. For example, isocurvature 
fluctuations in the baryons can be produced by inhomogeneous baryogene- 
sis (Turner, Cohen, & Kaplan 1989; Yokoyama & Suto 1991 and references 
therein), or by primordial magnetic fields (Kim, Olinto, & Rosner 1995). 
IHDM models were apparently first proposed by Peebles (1983) and inves- 
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tigated more recently by Sugiyama, Sasaki, & Tomita (1989). In light of 
the stringent constraints which the COBE results place on any model for 
large-scale structure, this is an opportune time to re-examine these models. 

In the next section, we discuss the initial conditions for the IHDM model, 
and calculate the cosmic microwave fluctuations and power spectra for isocur- 
vature HDM models with an initial power-law fluctuation spectrum, normal- 
izing to the COBE results. In section 3, we examine mixed isocurvature and 
adiabatic HDM models (MHDM), in which the former component dominates 
on small scales and the latter on large scales. In section 4 we discuss the data 
from the cosmic microwave background (CBR) and structure formation with 
which we will constrain our models. Finally, in section 5 we investigate a 
unique characteristic of isocurvature HDM models: such models naturally 
produce a bias between baryons and HDM even in the linear regime. In 
section 6 we present the results of our fits to the data. Our conclusions are 
summarized in section 7. We find that models with only isocurvature per- 
turbations and an initial power law spectrum are ruled out by the existence 
of damped Lya systems, but models with a mixture of adiabatic and isocur- 
vature perturbations can fit the data; in particular, such models with white 
noise baryon fluctuations give an acceptable fit. 

2 Power spectra and microwave fluctuations 

We assume an ft = 1 universe dominated by a single massive neutrino with 
m v = 93(1 — Hsj/i 2 , where h is the Hubble parameter in units of 100 km 
sec -1 Mpc -1 and Qb is the baryon fraction. (This differs from the models 
examined by Sugiyama et al. (1989), who assumed three massive neutrinos of 
equal mass). We take h = 0.5, consistent with universal age constraints, and 
assume Qb = 0.0125/i -2 , consistent with primordial nucleosynthesis (Walker 
et al. 1991). 

The initial isocurvature fluctuations in the baryons are assumed to have 
Gaussian statistics and initial power spectrum Ps(k). Because the fluctua- 
tions are isocurvature, the photons and neutrinos initially have compensating 
fluctuations which cancel the baryon fluctuations to give a small net pertur- 
bation (for a detailed analysis see Schaefer and de Laix 1995). For simplicity, 
we take Pb(A;) to have a power-law initial spectrum: 



P B (k) = Ak 
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where the amplitude A and power index n# are to be fit to the observations. 
We calculate the final processed linear power spectrum P(k) using the for- 
malism of Schaefer and de Laix. For each component, this involves solving 
the linearly perturbed Boltzmann equation given by 



Here, T is the gauge invariant perturbation in the phase space density /, 
where T is a function of the comoving momentum v and wave vector k. T 
is the derivative with respect to conformal time, /i is the cosine of the angle 
between k and v and q is the comoving energy. The variables A, V and II are 
moments of the phase space density representing the energy density, velocity 
and anisotropic pressure fluctuations respectively (see Schaefer and de Laix 
1995 for more details). 

Asymptotic limits are easy to calculate: on scales larger than the hori- 
zon at equal matter and radiation, the usual isocurvature suppression gives 
P(k) oc k nB+4 , while in the limit of large A;, the free-streaming of the neutrinos 
produces a power spectrum which asymptotes to P(k) oc k nB ~ 4 (Villumsen 
et al. 1991). What are reasonable values for n#? If the process which pro- 
duced the fluctuations in the baryons was uncorrelated on scales larger than 
the horizon when the perturbations were generated, then we expect on those 
larger scales a white-noise spectrum, for which n B = 0. This spectrum is the 
Gaussian limit of the seeded hot dark matter model examined by Villumsen 
et al. (1991), who found that the model produced acceptable results on small 
scales. The limit n B = —3 corresponds to the isocurvature Zel'dovich power 
spectrum (one of the two cases examined by Sugiyama et al. 1989, who also 
investigated n B = —1). We will leave n B a free parameter to be fixed with 
respect to the observations contained in section 4. 

First let us consider the baryon transfer function Tb (k) defined by Tb (k) = 
Ab /Abo where Abo and Ab are the initial and evolved perturbations respec- 
tively. Thus, P B (k) oc T B {kfk nB . We are left with the freedom to choose 
a convenient normalization for Tg(A;). It is easily verified that the growth 
factor for baryon fluctuations outside the horizon is the same regardless of 
scale, making the baryons a convenient reference point for defining transfer 
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functions. We thus can choose the normalization such that Tg(A;) = 1 in the 
limit of small k, while the transfer functions for the total perturbation and 
the HDM are defined in proportion to Tb(A;). This is illustrated in figure 1 
which shows a plot of the linear transfer functions for the baryon, HDM and 
total perturbations evolved linearly to the present. For k given in units of 
Mpc, the total linear transfer function is well fit by (better that 1%): 

T (k) = ^! 

' \l + a 1 k + a 2 k 2 + a 3 k 3 + a 4 k 4 + a 5 k 5 + a 6 k 6 + a 7 k 7 \ a& 

/ n . an . 012V 13 (-Oia\ (q , 

where the values for the coefficients a; are given in table 1 for redshifts 
z = 0, 3.2, 4 and 9. Similarly one may fit the baryon transfer function with 

T ,^ _ Qq (1 + ai$k + ai 6 k 2 ) 

B ' ' ~ 1 1 + «i k + a 2 k 2 + a 3 k 3 + a 4 k 4 + a 5 k 5 + a 6 k 6 + a 7 k 7 \ a8 

«11 «12V 13 /— «14 



+ «.»( 1 + f + #J ^{-T 1 )- (4 ' 

where the coefficients are given in table 2. 

With these two transfer functions, one may derive the results for the hot 
component using the formula 

T H (k) = ^T(k)-^T B (k), (5) 

where the Si's are ratios of the energy density of the component to the closure 
density It is important for the reader to realize that these transfer functions 
do not necessarily scale in any simple way with redshift, unlike, for example, 
those produced by adiabatic CDM fluctuations which scale as + in the 
matter dominated epoch. To illustrate this point, we have plotted in figure 2 
the transfer functions for the baryons, hot component and total respectively 
evaluated at redshifts of 1 + z = 1,5, 10 and 20. For the case of the total 
density perturbation, it is often adequate to just scale the transfer function 
by (1 + z)^ 1 for moderate redshifts. However, one can see clearly from figure 
2 that such is not the case for baryons. 

We would now like to turn to the effects IHDM will have on the cosmic mi- 
crowave background. Because isocurvature perturbations produce only small 
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fluctuations on super-horizon scales, the causes of the fluctuations in the cos- 
mic microwave background are distinctly different from those produced by 
adiabatic models. For the case of adiabatic perturbations on large angular 
scales, like those probed by COBE, microwave distortions are generated by 
the Sachs- Wolfe effect, i.e. by the variations in the gravitational potential 
induced by density fluctuations (see e.g. Peebles 1993). Conversely, isocur- 
vature perturbations produce only small total fluctuations so that it is the 
intrinsic fluctuation of the photon density which dominates, not the Sachs- 
Wolfe contribution. While the sources are different, the observed multipole 
moments, at least on large scales, have a similar £ dependence, where £ is the 
order of the multipole. However, the amplitude for AT/T given similarly- 
normalized power spectra will be different. One can see the similarity in 
the low order multipole moments by looking at the temperature fluctuations 
in Fourier space. For perturbations generated by the Sachs-Wolfe effect, 
the Fourier component of AT/T with wavenumber k is related to the total 
density fluctuation A by AT/T oc A/k 2 . As it turns out, for IHDM, the 
perturbation in the photon density is also related to the total density by 
A 7 oc AT/T oc A/k 2 . Since AT/T scales similarly with k for each case, 
both adiabatic and isocurvature models will produce similar multipole dis- 
tributions on large scales, i.e. small £. On scales smaller than the horizon 
at recombination, the story is different. Here the microphysical interactions 
between the photons and baryons prior to decoupling determine the shape 
of the multipole distribution. The location and amplitude of Doppler peaks 
depends significantly on the velocity at which the perturbation in the photon 
baryon fluid crosses the horizon (see e.g. Hu and Sugiyama 1995); adiabatic 
perturbations cross the horizon with a significant velocity perturbation while 
isocurvature perturbations do not. In figure 3, we show a plot of the first 
1500 coefficients to the multipole expansion, ct£, for IHDM models with power 
indices running from n# = —3 to 0, where all coefficients are plotted as the 
ratio aj£(£ + 1) /90ag. For comparison, we also show the results for adiabatic 
HDM. The n# = —3 IHDM model is comparable to the HDM plot as they 
are both produced by initially scale invariant spectra, and one can see the 
similarities on large scales (low numbered moments) between adiabatic HDM 
and IHDM with n# = —3. One can also see that the Doppler peaks for adi- 
abatic and isocurvature models are out of phase, a result anticipated by Hu 
and Sugiyama (1995). 



6 



3 Mixed adiabatic and isocurvature models 



For the sake of completeness, we would like also to consider models which add 
an adiabatic (i.e. curvature) component to the IHDM model. One can find 
motivation for this approach from inflation. There are many compelling ar- 
guments in favor of inflation (Guth 1981), and it naturally leads to adiabatic 
fluctuations with a k 1 power spectrum, regardless of the presence of isocur- 
vature fluctuations. In such a model, the isocurvature fluctuations must be 
produced either during or after inflation to avoid being erased by inflation; 
many models for producing isocurvature fluctuations during inflation have 
been proposed (Barrow & Turner 1981; Bond, Kolb & Silk 1982; Yoshimura 
1983; Linde 1985; Seckel & Turner 1985; Fukugita & Rubakov 1986; Kofman 
1986; Kofman & Linde 1987; Kofman & Pogosyan 1988; Turner, Cohen & 
Kaplan 1989). Gratsias et al. (1993) suggested the addition of adiabatic 
Zel'dovich fluctuations to the seeded hot dark matter model in order to en- 
hance the power on large scales. Lilje (1990) argued in the opposite direction, 
pointing out that adiabatic Zel'dovich HDM produces acceptable power on 
large scales but needs an additional fluctuation component on small scales. 

We assume, then, an adiabatic component with an initial Zel'dovich power 
spectrum, giving the present power spectrum 



Our numerical analysis has shown that for an adiabatic power spectrum 
defined by P A (k) = 1.2 x W A x kT\ (k) Mpc 3 that 

a 9 A = a 9 cobe when 

A\ = 1. For convenience we have redefined TA(k) so when A\ is equal to 
unity, we have a COBE normalized adiabatic spectrum. We add this to our 
present power-law isocurvature spectrum 



Recall that T(k) is the total present isocurvature transfer function given by 
equation (3). We will further assume that the processes that produce the 
adiabatic and isocurvature fluctuations are statistically uncorrelated which 
eliminates any cross terms in the power spectrum. That is, the total power 
spectrum may be written as the sum of the adiabatic and isocurvature com- 
ponents, i.e. P(k) = PA(k) + Pj{k). We are left with a model completely 
described by three free parameters, A\, Ai and n#. The constraints on these 



P A (k) = A.TK^k. 



(6) 




(7) 
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parameters are obtained from observational data which we shall discuss in 
the following section. 



4 Limits From COBE and Structure Forma- 
tion 



For both IHDM and MHDM models we have a number of free parameters 
which we would like to constrain with observations. Fortunately, data is avail- 
able on many different length scales to provide constraints. On the largest 
scales, near the horizon size, we have the microwave background fluctuations 
as measured by COBE. On cluster scales (8/z _1 Mpc) we have data from 
the QDOT IRAS galaxy survey, and on small scales (0.13/z -1 Mpc) we have 
observations of damped Lya systems. To constrain our models with this di- 
verse set of measurements we have elected to perform a \ 2 minimization on 
all data simultaneously to determine which parameters best fit the data as a 
whole. From this we can establish a likelihood that a particular model can 
explain observations on many different length scales. So let us now consider 
each of the data sets in more detail. 

We shall first turn to the COBE data. An analysis performed by Gorski et 
al. (1994) on the two year COBE data has shown for Sachs-Wolfe tempera- 
ture fluctuations generated by a power law density spectrum, i.e. P(k) oc k n , 
that the ninth coefficient of the multipole expansion is a constant for a broad 
range of power indices. We will take the Gorski et al. value of ag cobe = 
8.3 ± 0.7 /iK/T 70 with T 7 o = 2.726 yuK as the the COBE normalization data 
point for our fits. Using our numerical simulations, we have determined an 
analytic fit for ag as a function of n# for the IHDM models. Note that for this 
calculation we have defined the power spectrum by P(k) = A 2 k nB T 2 (k) Mpc 3 
where k is given in units of Mpc -1 . The fit for ag, valid for the range n# = — 3 
to 2 is given by 



where the coefficients are given in table 3. 

Gorski et al. have also constrained the value of the power index using the 
COBE data. For IHDM models, that is models with no adiabatic contribu- 
tion, this power index constraint translates into n# = —3.02 ± 0.36, so we 





8 



shall include this in our fits of IHDM models. The COBE observations prefer 
a Harrison-Zel'dovich spectrum which is reflected in the constraint on n#. 
For MHDM models, the large scale temperature fluctuations are produced 
by the adiabatic component which satisfies the power index constraint, so we 
will ignore it in this case. 

To make any further constraints, we must turn to structure observations. 
For scales larger than 8/i _1 Mpc, non-linear evolution is not significant, and 
any data that is available should correspond closely to the linear power spec- 
trum. Such data is available from an analysis performed on the QDOT IRAS 
survey (Efstathiou et al. 1990) by Feldman et al. (1994). They were able to 
estimate the power spectrum in the range of k = 0.02 to 0.2 h Mpc, well 
within the linear regime. However, since the data on which they based their 
observations are from redshift surveys, they are not equivalent to real space 
measurements, and thus corrections must be applied to relate the observed 
results to theoretical power spectra. Redshift surveys tend to show greater 
clustering than is actually present (Kaiser 1987), so a correction factor is 
introduced such that 



P(k) -> bj 



2 1 
1 + 3b~! + Wj 



P(k). (9) 



The parameter bj is a bias factor which represents the uncertainty to which 
the IRAS galaxies trace the underlying mass distribution. Measurement of 
the large scale peculiar velocity field from the QDOT data performed by 
Kaiser et al. (1991) gives the value of bj = 1.16 ± 0.21. We shall allow bj 
to be a free parameter to be best fit by our models. A second correction is 
required on small scales since the peculiar velocities of the galaxies tend to 
overshadow the clustering effect (Peacock 1992). This can be accounted for 
by multiplying by a damping factor 

pw - f ^f^w- do) 

where R v = 4 Ah' 1 Mpc. 

Where the QDOT analysis yields information on large scales, the ob- 
servation of damped Lya systems yield constraints on small scales. Recent 
observations have detected such systems out to redshifts of z ~ 4 (Storrie- 
Lombardi et al. 1995). They estimate the density fraction of the neutral hy- 
drogen in damped Lya systems, Cl gas , for a flat (Cl = 1) universe with h = .5 
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and find that Q gas = 2.2 ± 0.4 x 10~ 3 for z = 4 and Q gas = 3.3 ± 0.6 x 10~ 3 
for z = 3.2. Data for smaller redshifts show strong indications of evolution, 
which tends to reduce fl gas , so we will not consider this data. We will as- 
sume also that all the gas is neutral, ignoring the possibility that some might 
be ionized. Damped Lya systems are optically dense, allowing only surface 
ionization, so this is not an unreasonable assumption. 

The length scales corresponding to damped Lya have undergone highly 
non-linear evolution, but we can estimate the Lya mass function from the lin- 
ear power spectrum using the Press-Schechter technique (Press & Schechter 
1974; see e.g. Peebles 1993). In this formalism, the differential density frac- 
tion is given by 



where Y is the helium mass fraction, v c = A c /a and A c represents the critical 
collapse density. The density contrast a is the smoothed rms mass fluctuation 
given by 



where R is the smoothing length and P(k, z) is the power spectrum at redshift 
z. For damped Ly« systems, we take the minimum value for R necessary to 
produce an observed cloud to be RoLya = 0.13/i 1 Mpc (de Laix, Scherrer 
and Schaefer). This is our lower bound for integrating equation (11); we take 
an infinite upper bound. The value for A c is uncertain because it depends 
on the details of the collapse. If one uses a simple spherical collapse model, 
one would derive a value of A c = 1.69; however, fits to hydrodynamical 
simulations have produced results as low as A c = 1.33 (see e.g. Klypin et 
al 1994). We shall consider both possibilities when calculating our model fits 
to damped Lya systems. 



One of the unanswered questions in large-scale structure is the extent to 
which the baryon distribution traces the dark matter. This uncertainty is 
quantified in terms of the bias parameter 6, defined by 







5 Natural Bias 




(13) 
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where <r(R) and <Jb(R) are the smoothed rms density fluctuations for the to- 
tal and baryon perturbations respectively. The possibility that b ^ 1 is usu- 
ally attributed to poorly-understood processes which operate in the strongly 
nonlinear regime, particularly effects due to pressure which affect only the 
baryons. The IHDM model is unusual in that the distribution of baryons and 
HDM is naturally biased in the linear regime. This is due to the fact that 
when the neutrino free-streaming length is large, the neutrino perturbations 
on small scales are erased, while the baryon perturbations remain; we shall 
refer to this effect as "natural bias" . However, as the neutrino free-streaming 
length drops, the combined neutrino and baryon perturbations can begin 
growing again, and this unbiased growing component eventually swamps the 
initial pure baryon fluctuation. Hence, we expect b to be largest on small 
scales, and to be a decreasing function of time. This effect is obvious in the 
power spectrum derived by Sugiyama et al. (1989), who did not comment 
on it. 

We can see from figure 1 that natural biasing can be significant on the 
scales of damped Lya systems, ~ 0.1 hr x Mpc. It is important then to see 
how this natural biasing will affect limits derived from damped Lya clouds. 
One would like to know the ratio of ps/ p in a collapsed cloud which, because 
of natural biasing, should not in principle be Qb- We consider here only 
the natural bias in the linear density field, ignoring the non-linear effects. 
(Of course, this problem can only be completely resolved with hydrodynamic 
simulations, which are beyond the scope of this work.) With this restriction, 
we shall define the enhanced baryon fraction which we shall use in our Press- 
Schechter integral as 



This formula gives the ratio of baryon density to total in the linear regime 
for an object smoothed on a scale R with an overdensity A c . One may be 
a bit suspicious that we have used the natural bias b derived from the rms 
fluctuation when in fact we are only considering objects in highly over-dense 
regions, where the value of b from natural bias might be quite different. 
To test this we have numerically generated Gaussian density fields for the 
baryonic and total perturbations smoothed over various scales for which nat- 
ural biasing is significant. Our results indicate that, when considering only 
overdense regions with A > A c under the condition A c /cr > 1, the ratio of 




(14) 
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Ab /A is on average equal to b. The effect of natural bias then is to modify 
the Press-Schechter formula to give the following: 



•gas — 



1 + Y V 1 + A c 





— diner v r e 



(15) 



We use this to estimate our natural biased density in damped L a systems. 
The value of b is highly sensitive to n#. For example, with R = 0.1 hr x Mpc, 
we find that b = 4.5 for n# = 0, while b = 1.1 for n# = —3. 



We have proposed two classes of models and have presented the data with 
which we can constrain them, and we now present the results from our various 
fits. First let us consider the simpler of the two models, pure IHDM with 
a power-law power spectrum. Since we have no adiabatic contribution, we 
are left with three free parameters (A2, iib and bj). We find the best fit 
values by performing a x 2 fit to the data discussed in section 4, specifically 
COBE, the QDOT power spectrum and damped Lya systems. The best fit 
iib is ~ —2, but this value of n# is a poor fit to the Lya data. It produces 
so little power on small scales (asymptotically P{k) ~ k~ 6 ) that no IHDM 
model consistent with COBE and QDOT can produce a significant density 
of damped Lya systems at high redshift. Because of this lack of small-scale 
power, changing the value of A c or including natural bias effects does not 
alter the fit. 

Now let us add an adiabatic component to the isocurvature. Figure 4 
presents the best fit values for Ai and n with all four parameters (A\, A2, ris 
and bj) left free assuming no natural bias and a value for A c = 1.69. Also 
shown are the 68% and 95% confidence limits plotted in the A2 — n# plane 
with A\ and bj held fixed at their best fit values. The value of A\ is fixed 
by COBE since the large scale microwave fluctuations are dominated by the 
adiabatic component for reasonable models. The adiabatic component also 
dominates the the power spectrum on large scales, so that it fixes the value 
of bj as well. The amplitude and power index (A2 and %) of the isocurva- 
ture component are then determined by fitting the limits from damped Lya 
systems but are restricted by the small scale results from QDOT as well. 
Figure 5 shows similar results as the previous figure accept that the value for 
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A c is now 1.33. Figure 6 shows the same plot as figure 4 with A c = 1.69 but 
now the natural bias effect discussed in section 5 is included. Finally, Figure 
7 shows the results including the natural bias effect with A c = 1.33. 

If we accept the lower bound on A c = 1.33 as reasonable, we see that 
white noise (n# = 0) models are allowed at the 68% confidence limit for 
both the natural biased and unbiased cases. These models are of particular 
interest as they would arise from any process which is uncorrelated on scales 
larger than the horizon, and may be considered in some sense "natural". 

Since we have considered all of the data simultaneously in our \ 2 if 
is important to examine each observation independently. First consider the 
COBE constraint. The normalization for the adiabatic HDM component, 
which determines the value for ag, is well within the 2-cr limit determined 
from COBE for both the natural biased and unbiased fits. Now consider 
QDOT. The large number of data points from the QDOT data set assures 
that any acceptable model should fit it well in a \ 2 analysis, a fact that can 
be verified visually. Figure 8 shows a plot of the QDOT data along with 
the best fit IHDM model and the best fit MHDM model with n B = and 
natural bias. Also included for comparison is the cold+hot dark matter model 
(C+HDM) with 25% HDM and 75% CDM. It is normalized to COBE using 
the some bj as MHDM. The bias needed to fit the QDOT data is given by 
bj fx 0.66 for either the natural biased or unbiased fit. This is 2.4<r from the 
mean estimate of bj = 1.16 ± 0.21 (Kaiser et al. 1991). Finally, consider the 
damped Lya systems. We would like to see how the theoretical predictions for 
flgas stand up to the observed data. For the best fit, unbiased MHDM model, 
estimates of the gas fraction are ft gas = 1.0 x 10~ 3 ,2.4 x 10~ 3 for z = 4,3.2 
respectively. That is a 3<x difference for 2 = 4, the larger difference. With 
natural bias, the best fit gas fraction increases to fl gas = 1.3 x 10~ 3 , 2.8 x 10~ 3 
for 2 = 4, 3.2 respectively. This gives a 2.5<x difference at z = 4, a significant 
improvement over the unbiased model. For comparison, the C+HDM model 
with A c = 1.33 gives tt gas = 3.2 x 10~ 3 , 5.6 x 10~ 3 for 2 = 4, 3.2 respectively. 
C+HDM greatly overproduces damped Lya systems, but this is not as bad as 
it seems. We chose an infinite upper scale for our Press-Schechter integral, 
but a finite upper bound, determined from the rotation of spiral galaxies 
which damped Lya systems are believed to be the progenetors, may be more 
realistic (Mo & Miralda-Escude 1994; Kaufman & Chariot 1994). One can 
find ways to reduce the observed number of systems, but it is more difficult to 
increase that number. The C+HDM model is probably a better fit. Finally, 
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in figure 9 we present the multipole distribution for the best fit, natural bias 
model with a white noise (iib = 0) input spectrum. The dashed line shows 
the isocurvature contribution, the dot-dashed the adiabatic contribution and 
the solid line shows the sum. The results show that on COBE scales, this 
MHDM model in indistinguishable from standard adiabatic models; however, 
for moments £ > 100, there are significant deviations from pure adiabatic 
models which could be detected in future observations. 

7 Conclusions 

The results presented in section 6 give a clear picture of the types of models 
consistent with current observation. The most significant conclusion is that 
pure IHDM models are ruled out by the existence of high redshift damped 
Lya systems. If we ignore these objects completely, IHDM is still lies on tenu- 
ous ground with regards to COBE, as the best fit power index of n# = —2.23 
lies outside the 95% confidence limit of the COBE power index constraint. 
One can reasonably reject IHDM with a power law initial spectrum as a 
viable model. 

Conversely, MHDM should be considered a reasonable candidate. The 
largest strike against it is the unnaturally small value for bj required to fit 
QDOT. Measurements from galaxy redshift surveys like QDOT pose sig- 
nificant problems for ft = 1 models which do not contain a cosmological 
constant. It is difficult to reconcile COBE constraints on the power spec- 
trum amplitude with bj = 1 and QDOT. However, there is evidence that 
redshift surveys underestimate the large scale power. A recent analysis of 
the APM optical galaxy survey (Gaztanaga 1995) has yielded a value for 
<T 8h- 1 Mpc = 0-95 ± 0.07, significantly higher than those computed with bj = 1 
yielding (J 8h -i Mpc = 0.75 ± 0.098 (Peacock and Dodds 1994). For the best 
fit white noise, natural biased MHDM model, <T 8 /,-i Mpc = 1.1, while the 
C+HDM model gives cr 8 /j-i Mpc = 0.91. 

MHDM models are also at the fringe when it comes to the observation 
of damped Lya systems, especially at z = 4. However, there is again more 
breathing room than at first appears. Because of the statistical difficulties, 
we did not include the fact that the redshifts at which fl gas are observed 
are uncertain as well. The 1 — o uncertainties quoted from Storrie-Lombardi 
et al. (1995) for the redshifts are z = 41qJ and z = 3.2^'^. If the actual 
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redshifts are lower than the mean value, the damped Lya constraint would 
be significantly loosened. 

While adiabatic + isocurvature models do not provide the best fit to the 
data currently available, they are not in obvious contradiction to any of the 
data. The C+HDM model offers a superior fit to both QDOT and damped 
Lya system data, but its superiority is not compelling. Thus, MHDM models 
produce a natural alternative to C+HDM models for large scale structure. 
The MHDM models have the advantage of requiring only a single dark matter 
particle (which is known to exist), but they require an additional mechanism 
to generate isocurvature perturbations. 
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Figure Captions 

1 . The isocurvature hot dark matter transfer functions for the total (solid) , 
baryon (dashed) and hot dark matter (dot-dashed) perturbations evolved 
linearly to the present. The normalization is chosen such that Ts{k) = 

1 for small k. 

2. The evolution of the linear isocurvature hot dark matter transfer func- 
tions in redshift: hot dark matter (top), baryons (middle) and total 
(bottom). Plots are shown for 1 + ^ = 1,5,10 and 20. 

3. AT/T multipole ratio plots for isocurvature hot dark matter models 
iib = 0, —1, —2 and — 3 (solid curves) and for adiabatic hot dark 
matter (dashed curve). 

4. The Ai — n best fit along with the 68% and 95% confidence limits for 
the mixed (adiabatic + isocurvature) hot dark matter model including 
no natural bias effect; A c = 1.69, A\ = 0.80 and bj = 0.68. 

5. Same as 4 except A c = 1.33, A\ = 0.81 and bj = 0.67. 

6. Same as 4 except now the effects of natural bias are included; A c = 1.69, 
A 1 = 0.82 and bj = 0.65. 
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7. Same as 4 except now the effects of natural bias are included; A c = 1.33, 
A 1 = 0.82 and bj = 0.65. 

8. The QDOT power spectrum (Feldman et al. 1994) along with cold + 
hot dark matter (solid), isocurvature hot dark matter (dashed) and 
adiabatic + isocurvature hot dark matter (dot-dashed) power spectra. 

9. The AT/T multipole moments for the best fit adiabatic + isocurva- 
ture hot dark matter model (including natural bias) with a white noise 
(iib = 0) initial baryon power spectrum. The dashed curve shows the 
isocurvature contribution, the dot-dashed curve the adiabatic contri- 
bution and the solid curve the sum. 

Table Captions 

1. Fitting coefficients for the isocurvature hot dark matter total transfer 
function T(k) given at selected redshifts, where k is given in units of 
Mpc -1 and h = 0.5. Entries are written in the form a(b) = a X 10 b . 

2. Same as table 1 for the baryon transfer function Ts(k). 

3. Fitting coefficients for cig as a function of n. Entries are written in the 
form a(b) = a x 10 6 . 
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Table 1: 



Coefficient 


Z = 


z = 3.2 


z = 4 


z = 9 


(ii 


-1.658 


-2.284 


-1.601 


-2.051 


a 2 


9.223(3) 


9.390(3) 


9.212(3) 


9.225(3) 


a 3 


-1.774(5) 


-1.718(5) 


-1.656(5) 


-1.4812(5) 


0,4 


3.100(6) 


2.943(6) 


2.852(6) 


2.420(6) 


a 5 


-1.554(7) 


-1.269(7) 


-1.246(7) 


-7.215(6) 


a 6 


3.720(7) 


1.739(7) 


1.945(7) 


-9.235(6) 




2.336(8) 


3.122(8) 


2.960(8) 


3.778(8) 


a 8 


5.388(-l) 


5.322(-l) 


5.338(-l) 


5.296(-l) 


Clg 


2.383(5) 


5.673(4) 


4.764(4) 


2.382(4) 


a 10 


1.283(-1) 


8.888(-2) 


8.577(-2) 


7.664(-2) 


a n 


-5.998(-l) 


-7.272(-2) 


7.206(-4) 


3.045(-2) 


a\2 


3.888(-2) 


9.809(-3) 


4.437(-3) 


1.274(-3) 


«13 


3.051 


4.334 


6.375 


1.026(1) 


«14 


2.003(-l) 


2.303(-l) 


2.839(-l) 


3.325(-l) 
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Table 2: 



Coefficient 


z = 


z = 3.2 


z = 4 


z = 9 




3.975(1) 


5.321(1) 


4.567(1) 


6.954(1) 


0,2 


8.330(2) 


4.552(2) 


8.306(2) 


-3.230(2) 


a 3 


-1.487(4) 


-1.378(4) 


-2.411(4) 


1.636(4) 


0,4 


3.791(5) 


7.804(5) 


8.224(5) 


7.328(5) 


a 5 


-2.368(6) 


-5.943(6) 


-5.847(6) 


-6.750(6) 


a 6 


7.555(6) 


1.982(7) 


1.861(7) 


2.547(7) 


0,7 


-4.190(6) 


-1.109(7) 


-1.035(7) 


-1.453(7) 


Clg 


1.087 


1.053 


1.073 


9.697(-l) 


Clg 


1.002 


1.001 


1.001 


1.001 




1.914 


1.668 


1.635 


1.525 


a n 


8.899(-l) 


3.919(-1) 


3.696(-l) 


2.7212(-1) 


0-12 


2.339 


7.210(-1) 


6.024(-l) 


4.178(-1) 


0-13 


1.321 


1.365 


1.419 


1.254 


<2i4 


3.699(-l) 


3.224(-l) 


3.381(-1) 


2.800(-l) 




1.605(1) 


4.268(1) 


3.93(1) 


5.984(1) 


a 16 


2.718(5) 


6.977(4) 


5.767(4) 


3.042(4) 



Table 3: 



Coefficient 






-13.66 


Cl 


-1.885 


c 2 


.2698 


c 3 


-4.238(-2) 


c 4 


-1.537(-2) 


c 5 


4.856(-3) 


c 6 


1.356(-3) 
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